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Objective: Single nucleotide polymorphisms (SNPs) accumulated frequently in the mitochon- 
drial displacement loop (D-loop) in many cancers. We had identified cancer risk-associated 
SNPs in the D-loop of non-Hodgkin lymphoma (NHL) patients previously, in this study, we 
investigated the association of age at onset and D-loop SNPs in NHL patients. 
Materials and methods: The D-loop region of mtDNA was sequenced for 133 NHL patients 
recorded at the Fourth Hospital of Hebei Medical University. The Kaplan-Meier method was used 
to identify age at onset-associated SNPs in the D-loop of NHL patients. The Cox proportional 
hazards model was used to identify independent risk factors for age at onset. 
Results: The SNP sites of nucleotides 146C/T, 151T/C, 194T/C, 315C/C insert, 523Del/A, and 
525Del/C were identified for their association with age at onset, by the logrank test. In an overall 
multivariate analysis, allele 146 (relative risk, 0.403; 95% confidence interval [CI]: 0.182-0.895) 
(P = 0.026), allele 151 (relative risk, 0.378; 95% CI: 0.165-0.868) (P = 0.022), and allele 315 
(relative risk, 3.554; 95% CI: 1 .344-9.400) (P = 0.01 1 ) were identified as independent predictors 
for age at onset in NHL patients. 

Conclusion: SNPs in the D-loop can predict age at onset in NHL patients. Analysis of the 
D-loop SNPs can help identify NHL patient subgroups at high risk of early onset. 
Keywords: Kaplan-Meier; Cox proportional hazards model, D-loop, single nucleotide 
polymorphisms 

Introduction 

Non-Hodgkin lymphoma (NHL) is the fifth most frequent form of cancer among men 
and women, accounting for 2. 1 5% of all cancer-related deaths in the People's Republic 
of China over the period 2003-2007. 12 Although the etiology of NHL is poorly under- 
stood the immune-related conditions, family histology of hematopoietic malignancy, 
and body mass index have been shown to have an association with the cancer risk. 3,4 
Some studies have also shown that genetic factors are important to development of 
this cancer. 5-7 The polymorphism of genes in the oxidative stress pathway also was 
identified for their association with a risk of NHL. However, the true mechanism of 
this cancer remains unknown. 8-10 Although some studies have focused on genetics 
associated with the age at onset of NHL, they have yet to demonstrate a genetic basis 
of this disease. 11 

The human mitochondrial genome is 16 kilobase (kb) in length and is formed as a 
closed-circular duplex molecule. 12 Mitochondrial deoxyribonucleic acid (mtDNA) is 
believed to be more susceptible to DNA damage and so acquires mutations at a higher 
rate than nuclear DNA because of the high levels of reactive oxygen species (ROS) 
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generation, lack of protective histones, and limited DNA 
repair capacity in the mitochondria. 1314 

The mitochondrial displacement loop (D-loop) region is 
important for regulating both replication and expression of the 
mitochondrial genome because it contains the leading-strand 
origin of replication and the main promoter for transcription. 
In cancer patients, sequence changes accumulated extensively 
in this region. 15-18 Previously, we sequenced the D-loop 
region in a case-control study of NHL patients and identi- 
fied 140 single nucleotide polymorphisms (SNPs), including 
26 SNPs with frequency distribution of minor allele greater 
than 5%, and we also identified NHL risk-associated SNPs. 19 
In this study, we investigated the association of age at onset 
and SNPs of the D-loop, in NHL patients. 

Materials and methods 

Tissue specimens and DNA extraction 

Diagnosis of lymphoma was made according to the WHO 
Classification of Tumors. 20 Blood samples were collected 
from 133 NHL patients who received treatment in the Respi- 
ratory Department at the Fourth Hospital of Hebei University, 
Shijiazhuang, People's Republic of China, between 1990 
and 201 1. The genomic deoxyribonucleic acid (DNA) was 
immediately extracted using the Wizard® Genomic DNA 
Purification Kit (Promega Corporation, Fitchburg, WI, USA) 
and stored at -20°C. All procedures were supervised and 
approved by the Human Tissue Research Committee at the 
hospital. Informed consent was obtained from all participants 
before enrollment. 

Polymerase chain reaction (PCR) 
amplification and sequence analysis 

A 982-base pair (bp) product of D-Loop was amplified with 
the forward primer 5'-CCCCATGCTTACAAGCAAGT-3' 
(nucleotide 16190-16209) and reverse primer 5'-GCTTT- 
GAGGAGGTAAGCTAC-3' (nucleotide 602-583). PCR was 
run with a PCR Master Mix Kit (Promega Corporation) and 
purified prior to sequencing. Cycle sequencing was performed 
with the BigDye Terminator v3. 1 Cycle Sequencing Kit (Life 
Technologies, Carlsbad, CA, USA), and the products were 
then read with the ABI PRISM® 3 1 00 Genetic Analyzer (Life 
Technologies). Polymorphisms were confirmed by repeated 
analyses from both strands. 

Statistical analysis 

The age at onset curve of NHL patients was calculated, using 
the Kaplan-Meier method at each of the SNP sites, and 



compared with the logrank test. Multivariate analysis was 
performed with a Cox proportional hazards model. All of 
the statistical analysis was done with the SPSS Statistics 
for Windows, Version 17.0 software package (SPSS Inc, 
Chicago, IL, USA). A P-value of <0.05 was considered 
statistically significant. 

Results 

The age at onset distribution of the NHL patients was listed 
in Table 1. Those analyzed included: 25 patients at age less 
than 30 years old; 16 patients at age 30-40; 28 patients at age 
41-50; 30 patients at age 51-60; 20 patients at age 61-70, 
and 14 patients at age older than 70. Age at onset was ana- 
lyzed together with clinical characteristics, including gender 
and TNM (tumor node metastasis) classification, using the 
Kaplan-Meier method, and these were compared by the 
logrank test. TNM stage was associated with the age at onset, 
at statistically significant levels (Table 2 and Figure 1). 

Twenty-six SNPs with minor allele frequency higher 
than 5% were subsequently evaluated for their association 
with age at onset. The NHL patients were divided into two 
groups on the basis of their genotype at each SNP site, 
and the association of SNPs in the D-loop and age at onset 
was calculated with Kaplan-Meier method. The SNP sites 
of nucleotides 146C/T, 151T/C, 194T/C, 315C/C insert, 
523Del/A, and 525Del/C were identified for their association 
with age at onset by the logrank test (Table 3 and Figure 2); 
the SNP site of nucleotide 16261 was also identified for its 
association with age at onset but with borderline statistical 
significance (P = 0.083). We performed multivariate analysis 
for these predictors, including these SNPs and the clinical 
characteristics, with the Cox proportional hazards model. 
Allele 146 (relative risk, 0.403; 95% confidence interval 
[CI]: 0.182-0.895) (P = 0.026), allele 151 (relative risk, 
0.378; 95% CI: 0.165-0.868) (P = 0.022), and allele 315 
(relative risk, 3.554; 95% CI: 1 .344-9.400) (P= 0.01 1) were 



Table 1 


Onset age distribution in NHL patients 




Age 


Cases (N) 




Male 


Female 


<30 


20 


5 


30-40 


5 


1 1 


41-50 


16 


12 


51-60 


20 


10 


61-70 


10 


10 


>70 


10 


4 



Abbreviation: NHL, non-Hodgkin lymphoma. 
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Table 2 Clinical characteristics and their association with onset 
age in NHL patients 



Characteristics 


erases (in) 


P-value 


Sex 




0.443 


Male 


81 




Female 


52 




TNM classification 




0.031 


l-ll 


43 




111— IV 


90 





Abbreviations: NHL, non-Hodgkin lymphoma; TNM, tumor node metastasis. 
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Figure I Comparison of age at onset for NHL patients according to (A) sex and 
(B) TNM classification with Kaplan-Meier methods. 

Abbreviations: NHL, non-Hodgkin lymphoma; TNM, tumor node metastasis. 



Table 3 Polymorphic sites of the D-loop and their association 
with onset age in NHL patients 



Polymorphic sites 


Cases (N) 


P-value 


1 46 




0.044 


C 


13 




T 


1 20 




I5I 




0.006 


T 


7 




C 


1 26 




1 94 




0.082 


T 


7 




C 


1 26 




1 626 1 




0.083 


T 


9 




C 


1 24 




3I5 




0.024 


C 


6 




C insert" 


1 27 




523 




0.00 1 


Del 


43 




A 


90 




525 




0.000 


Del 


46 




C 


87 





Note: including C and CC insertion. 

Abbreviations: D-loop, displacement loop; NHL, non-Hodgkin lymphoma. 



identified as independent predictors for age at onset in the 
NHL patients (Table 4). 

Discussion 

Many SNPs in the D-loop region have been identified for 
their association with cancer risk and disease outcome in 
cancers. 21-27 The present study has extended those analyses 
to determine the relationship between age at onset and 
germline SNPs, in NHL patients. The alleles 146, 151, 
and 315 were identified for their association with age at 
onset, at statistically significant levels, by multivariate 
analysis. 

We have previously identified age at onset-associated 
SNPs of the D-loop, in patients with hepatocellular car- 
cinoma and esophageal squamous carcinoma. 28,29 In this 
study, we also identified age at onset-associated SNPs in 
NHL patients. 

All of the cancer risk-associated SNP sites are located 
in the hypervariable segment (HVII) region, mutational 
hotspots at which germline and tumor mtDNA mutations 
preferentially occur. 30 We have identified outcome-associated 
SNPs of these regions in other cancers. 21,27 The allele 315 
was identified as a cancer risk-associated SNP for NHL, and 
allele 151 was identified as a cancer risk-associated SNP for 
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Figure 2 Comparison of age at onset for NHL patients, according to the genotype of nucleotides 1 46, 1 5 1 , 1 94, 1 626 1 , 3 1 5, 523, and 525 in the mitochondrial D-loop, using 
the Kaplan-Meier method. 

Abbreviations: D-loop, displacement loop; NHL, non-Hodgkin lymphoma. 



non-small cell lung cancer in previous study. 19 - 31 Our data 
imply the important role of HVII region in modifying the 
onset age of cancer. 

Since the mitochondrial D-loop is important for replica- 
tion and expression of the mitochondrial genome, SNPs in 
this region might alter the function of electron transport 
chains by modifying mtDNA replication and transcription, 
which is responsible for the generation of ROS, and could 
contribute to nuclear genome damage as well as cancer initia- 
tion and promotion. 32 35 These ROS-mediated mechanisms 
may accelerate earlier onset of this disease. 

In conclusion, SNPs in the D-loop can predict age at 
onset in NHL patients, and analysis of D-loop SNPs can 
help to identify NHL patient subgroups at high risk of 
early onset. 



Table 4 Multivariate analysis of predictors associated with age at 
onset of NHL patients 



Factors 


Relative risk 


95% CI 


P-value 


16261 


0.714 


0.329-1.548 


0.393 


146 


0.403 


0.182-0.895 


0.026 


151 


0.378 


0.165-0.868 


0.022 


194 


0.515 


0.1 92- 1.381 


0.187 


315-316 


3.554 


1.344-9.400 


0.01 1 


523 


1.330 


0.392-4.515 


0.648 


525 


0.400 


0.121-1.322 


0.133 


TNM classification 


1.601 


0.927-2.764 


0.091 


Abbreviations: CI, 


confidence interval; NHL, 


non-Hodgkin lymph 


oma; TNM, 



tumor node metastasis. 
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